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Constraining the QCD phase diagram by tricritical lines 
at imaginary chemical potential 
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We present unambiguous evidence from lattice simulations of QCD with three degenerate quark 
species for two tricritical points in the (T,m) phase diagram at fixed imaginary fj,/T = in/ 3 mod 
27r/3, one in the light and one in the heavy mass regime. These represent the boundaries of the 
chiral and deconfinement critical lines continued to imaginary chemical potential, respectively. It is 
demonstrated that the shape of the deconfinement critical line for real chemical potentials is dictated 
by tricritical scaling and implies the weakening of the deconfinement transition with real chemical 
potential. The generalization to non-degenerate and light quark masses is discussed. 
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The QCD phase diagram is at present largely un- 
known. It describes which different forms of nuclear mat- 
ter exist for different choices of temperature and baryon 
density, and whether they are separated by phase tran- 
sitions. Its knowldedge is thus of great importance for 
current and future experimental programs in nuclear and 
heavy ion physics as well as astro-particle physics. Since 
QCD is strongly coupled on scales of a baryon mass and 
below, fully non-perturbative calculations are warranted. 
Unfortunately, Monte Carlo simulations of lattice QCD 
at non-vanishing baryon density are prohibited by the 
"sign" problem. To date only indirect methods are avail- 
able, introducing additional approximations which are 
justified for fi/T<l only One of these consists of 
simulating QCD at imaginary chemical potential fi = ifii 
with /Xj 6 M, for which there is no sign problem, and ana- 
lytically continuing the results to real chemical potential 
@, HI ■ While the Monte Carlo results contain the full in- 
formation about imaginary /i, analytic continuation via 
truncated polynomials fitted to the data introduces the 
approximation. 

In this letter we propose instead to study the phase di- 
agram of QCD at imaginary chemical potential in its own 
right. We shall demonstrate that there are intricate first 
order, triple, critical and tricritical structures, whose de- 
tails depend on the number of dynamical quark flavours 
Nf and their respective masses m / . These structures are 
bona fide properties of QCD and for this reason alone 
merit a detailed investigation. Moreover, we show that 
tricritical lines found at imaginary chemical potential, 
with their associated scaling behaviour, represent impor- 
tant constraints for the critical surfaces at real chemical 
potential. Finally, the phase diagram we investigate may 
serve as benchmark for studies within effective models 
(such as PNJL, sigma models, quark hadron models etc.), 
which can be easily extended to imaginary \x. 

Here we present a study of the (T, m) phase structure 
of QCD at fixed imaginary chemical potential (/i/T) c — 



i(2n + 1)ttT/3, n = 0, ±1, ±2, . . . for N f = 3 degenerate 
quark flavours. At those values of //, QCD undergoes a 
transition between adjacent Z(3) sectors. This is due to 
the exact symmetries of the partition function, 
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for complex fi The different Z(3)-sectors can be dis- 
tinguished by the Polyakov loop 



1 Nt 

L(x) = -Trn^o(x,r) = |L|e- 
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whose phase tp cycles through (ip) = n(27r/3),n = 
0, 1, 2, . . . as the different sectors are traversed. 

Hence, for imaginary /i = ifii there is a global Z(3) 
symmetry, even in the presence of finite mass quarks. Its 
spontaneous breaking implies transitions between neigh- 
bouring sectors for the mentioned critical values. These 
are first order phase transitions for high temperatures 
and analytic crossovers for low temperatures 0-0] , as 
shown schematically in Fig. Q] (left). For fi = iirT the 
order parameter to distinguish between the phases is 
the imaginary part of the Polyakov loop, Im(L). At 
high temperature there is a two-phase coexistence with 
fluctuations about (Im(£)) = ±\/3/2, on either side of 
the boundary. At low temperatures Im(L) fluctuates 
smoothly between those values. The situation is identical 
for the other critical values of /Ltj, up to Z(3) rotations. 

Away from a critical value of fii, there is a chi- 
ral/deconfinement transition line separating high and 
low temperature regions. This line represents the ana- 
lytic continuation of the chiral/deconfinemcnt transition 
at real \x. Its nature depends on the number of quark 
flavours and masses. Early evidence 0, H| is consistent 
with this line meeting the Z(3) transition at its endpoint 
between first order and crossover, and our present analy- 
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FIG. 1: Left: phase diagram for imaginary p. Vertical lines 
are first order transitions between different Z(3)-sectors, ar- 
rows indicate the phase of the Polyakov loop. The = 
chiral/deconfinement transition continues to imaginary chem- 
ical potential, its order depends on Nf and the quark masses. 
Right: phase diagram for Nf — 3 at fixed n — inT. Solid lines 
are lines of triple points ending in tricritical points, which are 
connected by a Z(2)-line. 



sis unambiguously confirms this. The nature of the end- 
point of the Z(3) transition line has already been inves- 
tigated for Nf — 4 6] and more recently for Nf = 2 

In this letter we study the nature of this junction at 
fixed \i = inT in Nf = 3 QCD as a function of quark 
mass. There are three possibilities, which we shall find 
to be all realised. For small masses the chiral transition is 
first order and branches off the Z(3)-transition line, ren- 
dering the meeting point of the three first order lines a 
triple point. For intermediate masses, the Z(3) transition 
ends in a second order endpoint with 3d Ising universal- 
ity, i.e. the chiral/deconfinement transition in its vicinity 
is a crossover. For large masses there is a first order de- 
confinement transition meeting the Z(3) transition again 
in a triple point. Hence, for fixed \x — inT, we obtain a 
(T, m) phase diagram as in Fig. Q] (right). The endpoints 
of the solid lines, which separate triple points from Ising 
points, correspond to tricritical points. 

To establish the phase diagram Fig. [T] (right) numer- 
ically, we work on lattices with temporal extent Nt = 4 
with standard staggered fermions at fixed fi/T — in, us- 
ing the RHMC algorithm and setting aside possible issues 
with taking a fractional power of the fermion determi- 
nant. For fixed Nt, temperature is tuned by varying the 
lattice gauge coupling /?. Hence, for a given bare quark 
mass am, we investigate the nature of the transition as 
a function of f3. To determine this, we analyse the finite 
size scaling of the Binder cumulant 



B A {X) = {{X-{X)f)/{{X-{X)ff 



(3) 



with X = Im(L) and (X) = 0. For li/T = in, ev- 
ery /3-value represents a point on the phase boundary 
and thus is pseudo-critical. In the thermodynamic limit, 
B^{j3) = 3,1.5,1.604,2 for crossover, first order triple 
point, 3d Ising and tricritical transitions, respectively. 
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FIG. 2: Finite size scaling of B4 for a small quark mass. 
On the right, the critical exponent was fixed to v — 1/3, 
corresponding to a first order transition. 



On finite L 3 volumes the steps between these values are 
smeared out to continuous functions whose gradients in- 
crease with volume. The critical coupling /3 C for the end- 
point is obtained as the intersection of curves from dif- 
ferent volumes. In the scaling region around a critical /3 C , 
B4 is a function of x — (/3 — (3 c )L 1 / l/ alone and can be 
expanded as 



Bi(/3, L) = Bi(p c , oo) + ax + bx 2 + 



(4) 



up to corrections to scaling, with the critical exponent v 
characterising the approach to the thermodynamic limit. 
The relevant values for us are v — 1/3, 0.63, 1/2 when /3 C 
corresponds to a first order, 3d Ising or tricritical transi- 
tion, respectively. 

For each quark mass, we simulated lattices of sizes 
L = 8, 10, 12 (20 in a few cases), at typically 8-14 different 
(3- values, calculated £?4(Im(L)) and filled in additional 
points by Ferrenberg-Swendsen reweig hting 0. Fig. m 
shows an example for a light quark mass am = 0.04. 
£>4 moves from large values (crossover) at low /3-values 
(i.e. low T) towards 1 (first order transition) at large 
/3-values (i.e. high T). In the neighbourhood of the in- 
tersection point, we then fit all curves simultaneously to 
Eq. Q, thus extracting /3 C , Bi{(3 c , oo), v, a, b. We observe 
that the value of the Binder cumulant at the intersec- 
tion can be far from the expected universal values in the 
thermodynamic limit. This is a common situation: large 
finite-size corrections are observed in simpler spin mod- 
els even when the transition is strongly first-order 
Moreover, in our case, logarithmic scaling corrections will 
occur near a tricritical point since d = 3 is the upper crit- 
ical dimension in this case Q. Fortunately, the critical 
exponent v, which determines the approach to the ther- 
modynamic limit, is less sensitive to finite-size corrections 
and in Fig. [2] consistent with v = 1/3, its value at a first 
order transition. A check is to fix v to one of the univer- 
sal values and see whether the curves collapse under the 
appropriate rescaling, Fig. [2] (right). Note that the crit- 
ical coupling determined from the intersection of the B± 
curves Fig. [5] (left) is consistent with that extracted from 
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FIG. 3: Left: Critical exponent v at \xjT — in. Right: Dis- 
tribution of Im(L) at the endpoint of the Z(3) transition. 



the peak of the specific heat or the chiral susceptibility. 

Proceeding in this way, we have investigated quark 
masses distributed over a large range, with results sum- 
marised in Fig. |3] (left). We find unambiguous evidence 
for a change from first order scaling to 3d Ising scal- 
ing, and back to first order scaling as the quark mass is 
made larger. Note that, in the infinite volume limit, the 
curve would be replaced by a non-analytic step function, 
whereas the smoothed-out rise and fall in Fig. [3] (left) 
corresponds to finite volume corrections. 

The results from the finite size scaling of B4 can be fur- 
ther sharpened by looking at the probability distribution 
of Im(L) at the critical couplings p3 c , corresponding to the 
crossing points. This is shown in Fig. [3] (right) for masses 
am = 0.05, 0.1, 0.2, 0.3 for L = 16. The lightest mass dis- 
plays a clear three-peak structure, indicating coexistence 
of three states at the coupling /3 C , which therefore cor- 
responds to a triple point. The same observation holds 
for heavy masses. For am = 0.1,0.2 the central peak is 
disappearing and for am = 0.3 we are left with the two 
peaks characteristic for the magnetic direction of 3d Ising 
universality. 

Hence, for small and large masses, we have unambigu- 
ous evidence that the boundary point between a first 
order Z(3) transition and a crossover at fi = ittT corre- 
sponds to a triple point. This implies that two additional 
first order lines branch off the Z(3)-transition line as in 
Fig. Q] (left), which are to be identified as the chiral (for 
light quarks) or deconfinement (for heavy quarks) transi- 
tion at imaginary chemical potential. This is expected on 
theoretical grounds: for m = or +00, these transitions 
are first-order for any chemical potential. The fact that 
the endpoint of the Z(3) transition line changes its na- 
ture from a triple point at low and high masses to second 
order for intermediate masses implies the existence of two 
tricritical points. Our current data on Nt = 4 put these 
between 0.07 < amtdci < 0.3 and 0.5 < amt r ic2 < 1-5. 

Since our Nt — 4 lattice is very coarse, a ~ 0.3 fm, 
an important question concerns cut-off effects. These 
strongly affect quark masses, and thus in particular the 
values of the tricritical quark masses where the changes 
from a triple point to a critical Ising point happen. How- 
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FIG. 4: Order of the transition as a function of quark masses. 
Left: quark hadron transition at /J, = 0. Right: the Z(3)- 
transition endpoint at fi/T = in/3. 



ever, universality implies that critical behaviour is insen- 
sitive to the cut-off, as long as the global symmetries of 
the theory are not changed. Our calculation is therefore 
sufficient to establish the qualitative picture Fig. Q] (right) 
for the continuum phase diagram at [i = i(2n + l)nT/3 
for three degenerate flavours. 

Let us now discuss how this critical structure is embed- 
ded in the parameter space with non-degenerate quark 
masses, Fig. [4] (right). The case Nf — 3 corresponds to 
the diagonal, with two tricritical points separating triple 
points from second order points. With non-degenerate 
quark masses these tricritical points will trace out tricrit- 
ical lines, ml rlc (m u ,d)- In the case of heavy quarks the 
situation is qualitatively the same for any Nf = 1,2,3 
T5L 11 1 - In the light quark regime, there is an interplay 
between the Z(3) and chiral symmetries and the situa- 
tion may be more complicated. The findings reported in 
[j| imply the existence of a finite tricritical light quark 
mass also for Nf — 2. It would then seem natural that 
the tricritical points for Nf = 2,3 are continuously con- 
nected by varying the strange quark mass, though this 
needs to be confirmed by explicit calculations. We stress 
that all critical structure indicated in Fig. 0] (right) can 
be determined reliably with standard Monte Carlo tech- 
niques, and continuum extrapolations are feasible with 
current resources. Knowledge of this phase diagram in 
the continuum should provide valuable benchmarks for 
the description of the QCD phase diagram by effective 
models. 

In order to establish the connection between imaginary 
and real chemical potential, let us briefly recall the situa- 
tion at /i = 0, Fig. |4] (left). The deconfinement transition 
in pure gauge theory is first order. In the presence of dy- 
namical quarks, it weakens with decreasing quark mass 
until it disappears along the deconfinement critical line 
with 3d Ising universality. The critical point for Nf = 1 
was determined in llj, and more recently for Nf = 1, 2, 3 
in a strong coupling expansion for a coarse N t = I lat- 
tice [13] ■ Similarly, the chiral transition for Nf = 2 + 1 is 
first order and weakens with increasing quark mass, until 
it disappears at a chiral critical line with 3d Ising uni- 
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FIG. 5: Critical line m c (jj, 2 ) in the 3-state Potts model P3] 
(left) and for QCD in a strong coupling expansion [T3| (right). 



versality 12, 13| . When a chemical potential is switched 
on, these critical lines sweep out critical surfaces which 
continue in the imaginary fi (or — /i 2 ) direction and join 
the tricritical lines at /i = iirT/3. We shall now illustrate 
this for the deconfmement critical surface. 

By universality, the properties of a second order tran- 
sition are the same in a model sharing the same global 
symmetry, and for the dcconfincmcnt transition this is 
the 3d three-state Potts model with Hamiltonian 



;x),0(x+i) 



-^[WW + AVtx)], (5) 



where 0(x) is a Z(3)-spin, and the couplings are identi- 
fied as h = exp - (M - //) /T, h' = exp - (M + fi) /T, while 
k increases with temperature. The qualitative change of 
the critical deconfmement line with chemical potential 
can be calculated in this model and one finds the first 
order region to shrink with real /j [lij ]. The same obser- 
vation is made for full QCD on a coarse lattice with a 
strong coupling expansion [l(J ■ See Fig. [SJ 

Both calculations show the continuation of the critical 
mass to negative fi 2 and a non-analyticity when joining 
the i?(3)-transition at \i = inT/3. However, in both 
cases it has not been fully realised that this junction is 
tricritical. Since chemical potential enters the partition 
function of the Potts model in the same way as in QCD, 
it features the same Z(3)-transitions. We can therefore 
directly check our QCD results in the heavy mass region 
against those in the Potts model at the same value /i = 
inT/3. For the latter, the Binder cumulant of the spin 



magnetisation was measured [1J]. We reanalysed those 
data fitting them to the scaling form, Eq. (|4]), and indeed 
find a change from first order behaviour, v = 0.33, at 
large values of M/T to second order 3d Ising, v = 0.63, 
implying again a tricritical point. 

Quite generally, a tricritical point represents the con- 
fluence of two ordinary critical points, i.e. in the heavy 
mass region the critical endpoints of the deconfmement 
transition at ji = i7rT/3(l ± e) merging into the Z(2>) 
cndpoint at fi = iirT/3. The deviation from the symme- 
try plane, ((fi/T) 2 + (it/3) 2 ), corresponds to an external 
field in a spin model, and the way a critical line leaves a 
tricritical point in an external field is again governed by 
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Fig. [5] shows that the data from [l4| and [Io( can both 
be fitted to this form. An excellent description of the 
data is achieved, reaching far into the real chemical po- 
tential region. Thus we conclude that for heavy quark 
masses, the form of the critical surface of the deconfme- 
ment transition is determined by tricritical scaling of the 
Z(3) transition at imaginary [i = iirT/3. 

It is clear that the chiral critical surface will likewise 
terminate on the chiral tricritical line at /i = iirT/3. Un- 
fortunately, for the chiral critical surface no suitable effec- 
tive model for finite density is available and we presently 
do not know whether tricritical scaling also shapes the 
chiral critical surface. This could be answered conclu- 
sively by extensive simulations at intermediate values of 
Hi. However, we do find m c (/i = iirT/3) > m c (0), which 
indicates a monotonous reduction of the critical quark 
mass as fj? is increased. This independently confirms the 
weakening of the first order chiral transition with real 
chemical potential observed previously in 0, [l]! EH ■ 
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